Off-diagonal generalization of the mixed state geometric phase 
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The concept of off-diagonal geometric phases for mixed quantal states in unitary evolution is de- 
veloped. We show that these phases arise from three basic ideas: (1) fulfillment of quantum parallel 
transport of a complete basis, (2) a concept of mixed state orthogonality adapted to unitary evolu- 
tion, and (3) a normalization condition. We provide a method for computing the off-diagonal mixed 
state phases to any order for unitarities that divide the parallel transported basis of Hilbert space 
into two parts: one part where each basis vector undergoes cyclic evolution and one part where all 
basis vectors are permuted among each other. We also demonstrate a purification based experi- 
mental procedure for the two lowest order mixed state phases and consider a physical scenario for 
a full characterization of the qubit mixed state geometric phases in terms of polarization-entangled 
photon pairs. An alternative second order off-diagonal mixed state geometric phase, which can be 
tested in single-particle experiments, is proposed. 

PACS numbers: PACS number(s): 03.65.Vf, 42.50.Dv 



I. INTRODUCTION 



The concept of geometric phase, anticipated by Pan- 
charatnam 1] in his study of interference of classical light 
in distinct states of polarization and developed by Berry 
for cyclic adiabatic quantal evolution, has beeii_^ener- 
alized in several steps. Aharonov and Anandan 3] con- 
sidered the cyclic nonadiabatic case and pointed out that 
the geometric phase is due to the curvature of the space 
of pure quantal states. Samuel and Bhandari Q pro- 
vided a general setting for the geometric phase so as to 
cover noncyclic evolution and sequential projection mea- 
surements. Another line of development has been the 
extension of the geometric phase to the mixed state case. 
Uhlmann 5] was probably first to address this issue in the 
mathematical context of purification and Sjoqvist et al. 
discovered a mixed state geometric phase for nonde- 
generate density operators in noncyclic unitary evolution 
in interferometry. The parallel transport conditions in 
Refs. 0,0] have been shown to lead to generically differ- 
ent interference effects and the mixed state concept in 
Ref. 6J has been extended to the case of completely pos- 
itive maps as well as to degenerate density operators 
0. An experimental test of the proposed mixed state 
geometric phase in Ref. in the qubit case has been 
reported recently [Tol | , using nuclear magnetic resonance 
technique. 

The noncyclic geometric phases in Refs. 0,0] become 
undefined when the interfering states are orthogonal and 
the interference visibility vanishes. This leads to an inter- 
esting nodal point structure in the experimental parame- 
ter space that could be monitored in a history-dependent 
manner 0, ■ The existence of nodal points also led 
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Manini and Pistolesi [TjI to introduce the off-diagonal 
geometric phases for pure states in adiabatic evolution. 
These phases may carry interference information at the 
nodal points of the standard geometric phase. The adi- 
abaticity assumption of Ref. 11^ was subsequently re- 
moved by Mukunda et al. [iJI ^^^d the second order 
off-diagonal pure state geometric phase was verified by 
Hasegawa et al. 0| for neutron spin. More recently, 
the off-diagonal generalization of the mixed state phase 
in Ref. 6] for parallel transporting unitarities has been 
proposed and extended to general unitarities by 
the present authors. 

In this paper, we wish to elaborate further on the con- 
cept in Ref. [131 ■ Our first concern is to develop sys- 
tematically the off-diagonal mixed state geometric phase 
in unitary evolution, filling in some important concep- 
tual gaps of Ref. [13 • The theory of parallel transport 
of a complete basis is considered. Under such a paral- 
lel transport, it is argued that the mixed state geometric 
phase follows naturally from a concept of orthogonality 
between density operators adapted to unitary evolution 
and a certain normalization condition. Our second con- 
cern is to provide a method for computing mixed state 
geometric phase factors to any order for the important 
class of unitarities under which the parallel transported 
basis of Hilbert space is divided into two parts: one part 
where each basis vector undergoes cyclic evolution and 
one part where all basis vectors are permuted among each 
other. The experimental realization of the two lowest or- 
der mixed state phases is discussed, using the idea of 
purification of the involved mixed states by attaching an 
ancilla system entangled to the considered system. We 
further develop the two-photon experiment proposed in 
Ref. ^3 so as to include also the lowest order mixed 
state phase. Finally, as an attempt to avoid the apparent 
need for an ancilla system, we consider an alternative off- 
diagonal mixed state geometric phase conceptually based 
upon the neutron experiment in Refs. |l5lll6j|. 
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II. BASIC IDEAS 
A. Quantum parallel transport 

Consider a continuous one-parameter family {U {s),s G 
[so, Si]|i7(so) = /} of unitarities that map any initial 
complete orthonormal basis {|V'fc)} of a Hilbert space H 
of dimension A'' to a continuous set of complete orthonor- 
mal bases {ji/'feCs)), s e [so,si]} of the same H. In the 
context of off-diagonal phases, it proves useful to consider 
parallel transport of such a complete set of orthonormal 
pure state vectors. 

Parallel transport in terms of the initial basis {|'0/c)} 
may be formulated as follows. Let J{s) be the Hamilto- 
nian operator in the Heisenberg representation. For the 
corresponding evolution operator U{s) we have {h ~ 1 
from now on) 



iU{s) ^ U{s)J{s), 
the formal solution of which reads 



U(s) 



V^'^ exp 



i / J{s')ds' 



(1) 



(2) 



Here, V^^ is inverse path ordering, i.e., ordering increas- 
ing s from left to right in each term of series expan- 
sion of the evolution operator. U{s) is said to parallel 
transport the initial basis {IV^fc)} if the local accumula- 
tion of phase along the unitary path vanishes for each 
\ipk), which amounts to (?/;fc|[/'l'(s)[/(s)|V'fe) = 0, i.e., 



(^fclJ(s)l^fc) =0, Vfc. 



(3) 



Thus, the one-parameter family of Hermitian generators 
J(s) has to be off-diagonal in the parallel transported 
initial basis and therefore traceless in any basis. In other 
words, U € SU(N) is a necessary condition for U being 
parallel transporting a complete basis. The converse does 
not hold: there are SU(N) transformations that contain 
generators that have nonvanishing diagonal elements in 
the {IV'fc)} basis. 

We may equally well formulate the parallel transport 
conditions in terms of the instantaneous basis {\'ij}k{s))} . 

Here, we have (?/'fe(s)|[7(s)f/^(s)|V'/c(s)) = 0, which en- 
tails 



{Ms)\H{s)\Ms))=Q, Vfc, 



(4) 



where H{s) — U{s)J{s)W{s) is the Hamiltonian opera- 
tor in the Schrodinger picture. Thus, H{s) has to be off- 
diagonal in the instantaneous parallel transported basis, 
which is consistent with U G SU(N). 

Any parallel transporting unitarity is denoted by f/" in 
the following. Moreover, an initial (instantaneous) non- 
degenerate [13 density operator whose eigenvectors co- 
incide with the basis ({l^fc (*))}) is said to be the 
parallel transported by C/" fulfilling Eqs. ijSJl and Q. 

For dimTY — N, any J{s) fulfilling the parallel trans- 
port conditions may be written in terms of N'^ — N 



linearly independent off-diagonal generators in the ba- 
sis {|?Afc)}. As an example, consider the qubit case 
= 2. Here, we expect J(s) be dependent upon two 
of the Pauli operators. Let {IV'i)) |'02)} be a parallel 
transported basis. Then J(s) = a{s)ax + h{s)(jy with 
(Tx = |V'l>(V'2| + |V'2>(V'l|, (Ty = -i|V'i)(V'2| +i|V'2)(V'i|, 
and a(s), h{s) being scalar functions of s. Any density op- 
erator p of the form p = Ail^Ai) (V'll +'^2|V'2) (V'2|, Ai 7^ A2 
is parallel transported by the corresponding unitarity. 
The same holds true for H{s) by making the replacement 

{\^i)M2)}^{\Ms)).ms))}. 



B. Orthogonality 

Two pure state vectors are orthogonal if their scalar 
product vanishes. On the other hand, any useful scalar 
product between density operators does not have this 
simple property and the concept of orthogonality be- 
comes less straightforward in the mixed state case. In- 
stead, for a given density operator p, one may take an- 
other density operator p' to be orthogonal to p if it yields 
minimum of the Hilbert-Schmidt product Tr[pp'] or the 

Bures fidelity ^ Tb{p,p'\ = [Tr^^p'^J^ the lat- 
ter being a worst case measure of distinguishability be- 
tween p and p' 21]. Here, we take a third approach to 
the concept of orthogonality adapted to unitarily con- 
nected density operators. The idea is to say that p and 
p' = UpU^ are orthogonal whenever they cannot interfere 
in the sense of Ref. . 

To develop this idea in detail, let us first suppose 
and |(p) are Hilbert space representatives of two arbitrary 
pure quantal states ip and (p, and assume further that {-ip) 
is exposed to the variable U(l) shift e'-*^. The resulting 
interference pattern obtained in superposition is deter- 
mined by the intensity profile p3 | 



e'X|^)-f|(^) -2 + 2|(V^|(/,)|cos[x-arg(V'|(^)] (5) 



that oscillates as a function of x- The key point here 
is to note that ^ and ip are orthogonal if and only if X 
is independent of x so that the interference oscillations 
disappear. 

This feature translates naturally to the mixed state 
case. Consider a pair of isospectral nondegenerate den- 
sity operators 



^^Xk\ipk){ipk\, P.^ = ^ Afe|(^fe)(<pfe 



(6) 



where each \ipk) — U\4>k) for some unitarity U. Each 
such orthonormal pure state component of the density 
operator contributes to the interference according to Eq. 
©. Thus, the total intensity profile becomes 

X (X Y.Xk\e'^\iJk) + \v^k)f 

k 

= 2 + 2^AA;|(V'fe|(^fe)|cos [x-arg(V'fc|(pfc)], (7) 
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where we have used that the A's sum up to unity. Fol- 
lowing the above pure state case, we say that and 
are orthogonal if and only if I is independent of x for a-U 
Hilbert space representatives {IV'fc)} and {|</3fc)} of the 
eigenstates of p^ and p,^, respectively. It follows that p^ 
and Pip are orthogonal if and only if {ipk\fk) = 0, Vfc. 

For an N dimensional Hilbert space Ti, we may gen- 
erate a set of N mutually orthogonal density opera- 
tors as follows. Assume pilfAfe) = A^IV'/c) is nondegen- 
erate and introduce a unitary operator Ug such that 

IV'n) = {Ug)^ n — 1, . . . , N . Thus, we may write 

Ug = |^l)(^jv| + IMi^N^ll + ■■■ |V'2)(V'l| (8) 

and it follows that 

Pn = {Ugyp,{Uiy, n = l,...,N (9) 

is a set of mutually orthogonal density operators. Ex- 
plicitly, this entails that 

Pi = Ai|?/'l)(V'l| + A2|'(/'2)(V'2| + ■ • • + AAr|V'Ar)(V'ivL 

P2 = Ai|V'2)(V'2|+A2|V'3)(^3l+--- + AAr|V'l)(V'l|, 
• ■ * ? 

PN = Ai|V'Ar)(^/'Ar| -I- A2|V'l)(V'l| + ■ ■ • 

+XnHn-i){^n-i\- (10) 

Notice here that different sets of mutually orthogonal 
mixed states may be generated by permuting the ■i/;„'s 
in Ug. 

C. Consistency and normalization 

The final step towards the concept of off-diagonal 
mixed state geometric phase, is to determine how the 
mutually orthogonal density operators should appear in 
the trace. This may be resolved as follows. 

We first notice that the Manini-Pistolesi expression 
[l3j| may be written in terms of pure state projectors 
Pjk = \i'jk){i'jk \ as 

7li^,P,,...P,, ^ $[Tr(C/llp,,C/llp,, . . . C/IIP,,)], (11) 

where $[z] = z/\z\ for any complex number z. We pro- 
pose to replace each of these projectors with F^^\pji^), 
where, for reason of permutation symmetry of the in- 
dexes ji , j2 , . . . , ji , the form of the function P*^'-* may only 
depend on Z. To assure consistency with Ref. we fur- 
ther require that P^-'^Pifc) ~^ Pjk pure state limit. 
We take the simplest nontrivial choice fulfilling this re- 
quirement, which is F^'-\pj^) — (F-^ [23, V = p(0 and 
q = q{l) integers p^ . 

Next, from (Pfe) = Pfe, we obtain the normalization 
condition 

Tr(C^gPfeC/gP(fe+l) modW • ■ ■ UlP{k+l) modw) 

= Tr((C/t)'Pfe) = 6iN, Vfc e [1, N], (12) 



where we have used C/J defined in Eq. lO, P{k+n) mod n = 

(C/g)""^Pfc(C/])""^ and (C/t)^ = /. We propose to de- 
mand that this normalization structure is preserved in 
the mixed state case. After the replacement Pj^ (^j^ ■• 
we similarly have 

^^y-' gPk '-^gP{k+l)modN ■ ■ ■ ^ gt^(k+l) mod N ) 

= Tr((C/t)Vr'') = 5.^Tr(pf ^), Vfc G [l,iV], 

(13) 

where we have used that (UgpUlY^"^ = UgpP/'^Wg. Thus, 
only p{N) — 1 and q{N) = N assures the desired kind of 
normalization in the mixed state case. Since p and q are 
functions of I only, it follows that p — 1 and q — l. 

This choice may also be understood from the following 
simple convergence arguments in the N ^ oo case. p'P^i 
typically involves factors of the form A'^/'', < A < 1. If 
Ip/q < 1 {Ip/q > 1) then the trace diverges (goes to zero) 
when I — !■ oo. Thus, only for Ip/q = 1 the N —> oo limit 
is finite and well defined. 



III. OFF-DIAGONAL MIXED STATE 
GEOMETRIC PHASE 

We are now ready to state our main result. The 
off-diagonal mixed state phase for an ordered set of 
I < N mutually orthogonal nondegenerate density op- 
erators pjf_, k = parallel transported by [/" is 
naturally given by 

^1'L,...P,,^'^MU\\^U\\^...U^\^)]. (14) 

This is manifestly gauge invariant and independent of 
cyclic permutations of the indexes ji , j2 , ■ ■ • , ji ■ By con- 
struction it reduces to Eq. Hll() in the limit of pure states. 
The mixed state geometric phase factor 

7« =$[Tr(C/llp,J] (15) 

proposed in Ref. 6] may be seen as a natural consequence 
of this general framework if we put 1 = 1. In Sec. V we 
propose experimental realization of this first (1 = 1) and 
second order {I = 2) phases, the latter being defined by 

in polarization-entangled two-photon intcrferometry. 

IV. COMPUTATION OF OFF-DIAGONAL 
MIXED STATE PHASES 

In the qubit case N = 2, consider the unitarity 

c/ii = c/|i I V'l) (^1 1 + C/I2 1 Vi)(V'2 1 + c/|i 1^-2) (^1 1 

+ C/|2l^2)(^2| (17) 
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that parallel transport some orthonormal basis 
{|^i>, IV'2)}. The matrix elements of C/H fulfill uf^ 
= (C/I2)* = rie-'^y^ and uf^U^^^ = -detf/H + uf^U^^^ = 
-l + r?2 as C/ll e SU(2). Here, t] = | (tAi If/" |^i> | is the 
pure state visibility and Q is the solid angle enclosed by 
the path traced out by the basis vectors {I"!/"!), IV'2)} and 
the shortest geodesic connecting its end points on the 
Bloch sphere. 

Now, in Eq. (|17|l parallel transports the mu- 
tually orthogonal density operators pi = Xi\ipi){Tpi\ + 
A2|V'2)('02| and p2 = Ai|V'2)(V'2| + X2\ipi){ipi\, for which 
we obtain 

Tr(C/llpi) = Tr(C/llp2)* 

Tr(i7llVprf/"VPJ) = -l + ??' + 2r7'v/A^cosl) 

= -1 + 7?^ +7?V-^b[pi,/02] cos 17, 

(18) 

where we have used the Bures fidelity J^b[pi,P2] ~ 

[Tr^^/prp2VPi] = 4A1A2. Notice that Tb[pi,P2] = 
for pure states and J-b[pit P2] = 1 in the maximally 
mixed state case. 

In the nondegenerate mixed state case Ai 7^ A2, the 
I = 1 phases are indeterminate only for 77 = 0, for which 
the I = 2 phase is well-defined since Tt(U^^ yfpiU^^ y/])2) = 
— 1. In the degenerate case Ai = A2, the density opera- 
tors pi and P2 become identical and spherically symmet- 
ric, so that no specific basis is singled out by the parallel 
transport condition and the mixed state geometric phase 
factors 7'"'^' and 7'-^-' therefore become undefined. Still, 
there is a unique notion of relative phase in this case with 
additional nodal points, as discussed in Ref. T^]. For a 
generic U — e"*'^"'"', cr = {ax, <^y,o-z) being the standard 
Pauli operators and |n| = 1, we obtain for Z = 1 nodal 
points at Tr(t/pi) = Tr(C/pi) = cos (5 = at which 5 

we have Tt{U ,/piU ^/p2) — cos 25 ~ —1. This shows 
that the I — 1 and / — 2 phases never become indetermi- 
nate simultaneously and thus provide a complete phase 
characterization of the qubit case. 

The off-diagonal mixed state geometric phase in the 
qubit case has a nontrivial nodal structure that arises 
due to the nonvanishing Bures fidelity. This can be seen 
by putting the left-hand side of Eq. (|18|l to zero and 
solving for ry^ yielding 

if = {1 + ^Tb[pi,P2] cosn)'\ (19) 

which has solutions at r; < 1 for Tb[pi, P2\cosVL > 0. 
Thus, the off-diagonal mixed state geometric phase fac- 
tor may change sign across the nodal surfaces in the pa- 
rameter space {J^b[pi, P2],Vi^) defined by the solutions 
of Eq. 1(1^ . as shown in Fig. ^ Thus, the corresponding 
off-diagonal mixed state geometric phase can take both 
values and tt, contrary to the corresponding pure state 
phase, which can only be tt. 
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FIG. 1: Nodal surfaces of the off-diagonal mixed state geo- 
metric phase for a qubit with the solid angle fl in steradi- 

ans. For Bures fidelity — |^Tr^yprp2ypi j > (mixed 
states), there are nodes also for paths with pure state visibility 

I II P 

rj = {?/>i|f/" 7^ 1 at various solid angles. 

To generalize the discussion we proceed to arbitrary 
Hilbert space dimensions N and provide a method for 
computing mixed state geometric phases to any order 
I < N for unitarities under which the parallel transported 
eigenbasis . . . , IV'at)} of the mutually orthogonal 

/o's is divided into two parts: one part where each basis 
vector undergoes cyclic evolution and one part where all 
basis vectors are permuted among each other. By appro- 
priate labeling of the eigenvectors, such unitarities can 
always be decomposed into the direct sum 

[/II = ujl (20) 

where ul permutes \i>m) ■ . . ^ IV'2) 

and m|| is diagonal in the remaining N — m cyclic eigen- 
vectors. These terms do not mix so that one may write 

lr(C/ll^...[/ll^) =7'(^)...,^^ (21) 

where P^'^^...p^., = Tr(u|^/p7r • • • 4^/ft7) and P^'^ -p., = 

Turning first to the contribution from the diagonal part 
of [/II , we have 

N 

E {uLy<A^r:^. (22) 

k—rn-\-l 

with ulf. the matrix elements of in the eigenbasis of 
the p's. As the density operators are nondegenerate, it 
follows that all Afc^ are different in each term on the right- 
hand side of Eq. H22|) and 'Dp^_^,,,pj^ must vanish if / > 
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rank of the p's. Notice here that arg Ulf. is the standard 
cychc geometric phase for the pure state "0^. 

Considering the contribution from Up we can estabUsh 
the following, li I = K x m, K integer < N/m and 
m > 2, then 



p(0 



.p, = [(-ir-Metu|]''/«...,^/Ai,...,A^),(23) 



where the /'^'^'s can be written as a sum of m terms. For 
other I, the P^'^'s vanish as there is K x m steps needed 
to connect ^pj^ and ^J/pj7 with u|. 



In the extreme case where all N eigenvectors are per- 
muted, only 7^p^^..p^^ may be nonvanishing. Here, K = 1 

1 since f/ll € SU(N). It follows 



and det Up — det Jj" = 
that 



Ph---Pn 



(-1) 



N- 



lr(N) 



(Ai 



(24) 



where each /'•^-' is determined by the sequence of p's. 
These /'s have some interesting properties. First, it can 
be seen that 



J Pi 



(N) 



1, 



(25) 



as a consequence of the normalization condition de- 
scribed in Sec. II. C. Thus, there exist at least one well- 
defined off-diagonal mixed state phase for [/H = uj,, in- 
dependent of the rank of the p's. Secondly, we have 

4^)..p,„ >0, Vji,...,jAr, (26) 

as each /^^'' always can be written as a sum of positive 
functions of the A's. This implies that the off-diagonal 
mixed state phases for such unitarities are completely 
determined by the dimension of the Hilbert space H. In- 
deed, for sequences where /^^^ 7^ we have 
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(JV) 



= —1, if dim(7i) even, 
-1, if dim(H) odd. 



(27) 



Let us now turn our attention to partial permutations 
characterized by m ^ TV, where we can use the follow- 



ing algorithm to determine /, 



(0 



As in the m — N 



case, each /p'| . 



iP3i---Pn 

■ ■Pii (Ai, • ■ • , Aat) decomposes into a sum of 
terms determined by the sequence of p's. Explicitly, we 
can write /p'|...pjj (Ai, . . . , \n) as a sum of m terms 



4:!...p,(Ai,...,A^) 



i=i 



A,(Ai 



(28) 



where Aq is of the form y A"^^ ■ ■ ■ -^Ti with integers 
ranging from to /. Since we are not interested in the 
phase contributions from uj, that are already included 
in the factor [(—1)™^-'^ detup] , we replace Up with the 
operator 



IV'mXV'll + l'0l)(V'2| + • ■ ■ \lpm~l){lpr, 



(29) 



being unitary on the permuted subspace. Thereafter, we 
compute /^p^...p^^(Ai,...,AAr) - Tr{W3^...W^^) with 

'[yjk ^ Ug"^"^ ^pji^ . Applying ordinary matrix multipli- 
cation rules and noting that only one entry in each row 
and column of W^'' is nonvanishing, Ai is in index nota- 
tion given by 



mod m (i+l) mod m, (2+2) mod rti 
(2+/ — 1) mod m,2. 



(30) 



Each eigenvalue Afe of pi appears exactly once in each 
W^^ , so our aim is to describe the correspondence be- 
tween the components W^\, and Afe by appropriate index 
transformations. 

To revert the off-diagonal matrices W^^ to diagonal 
form we apply (J/g™'')^ to each VF^'=, thus W'^'' — > 

{\/Pjk) — (t^g™"*)^!^"""- This transforms the indexes ac- 
cording to 



X 

y 



x' = {x + I) mod m, 



y 



leading to 



(31) 
(32) 



' (a; + l) mod m.y 

in terms of components. Consequently we obtain 

— (\//'^)(i+l)modm,(''+l)modm ■ ■ ■ 

^-1) mod m.(i-\-l—l) mod m 



' — 1) mod 771,(2+^—1) mod m ' 

(33) 

The transformation back to the "unpermuted" pi can 
be achieved by applying Ug described in Eq. ac- 
complishing |i/)„) — {Ug)"' {ipi), n = l,...,iV, thus 

^ ^ ^/P^ = iUly'-'^iUgV'-'- Since jk - 1 steps 
are needed to convert -^Pj^. to ^J/pT we have to carry out 
the index transformation 



x' ^ {x~ {jk - 1)) modiV 



(34) 



so that {^/pK)a:.x ~* iy^)x'.x'- ^ince pi is diagonal 
with eigenvalues Afe in ascending order in fc, the index x' 
denotes the wanted eigenvalue A^;'. 

This algorithm traces the locations of the eigenvalues 
Afe from back to ^J/pT when the unitary transforma- 
tions along the path starting from IF" are applied: 



(35) 



This analysis makes it possible to calculate the factor 
/pjj...Pj, (Ai, . . . , Aat) more efhciently than performing a 
multiplication of the I matrices involved. 
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Let us revisit the qubit {N = 2) case using the above 
general theory. If to = 0, both 7^^-' and jp^^ exist. More- 
over, we have 

T^^^L = ^^[KY + {uLy], (36) 

which is consistent with Eq. (|18|l for r] — 1. In the 
permutation case to = 2, we may use (—1)^^^ dot t/" = 

— 1 and, from Eq. if^ . /p^pa = 1 fo^' N ^ 2, and obtain 

Hfp2 = -1' (37) 

in agreement with Eq. H18|l for 77 = 0. 

As a further iUustration, let us work out the = 3 
case in detail. For to = 0, all the 7(^)'s are well-defined. 
The dependence upon the rank of the density operator is 
visible for higher I, namely 

^(3) ^ -r)(3) 

P1P2P3 P1P3P2 

(38) 

(2) (2) 

with Vp^pa and Vp^p^ obtained by permutations of the 
A's. In the to = 2 case, iV'i) IV'2) \ipi) while lips) 
undergoes cyclic evolution. Explicitly we have 

-^'^Lp. = T^fL,. = i/A^(t/|3)^ (39) 

(2) 

■PpiP2 can be calculated via the algorithm above as fol- 
lows. Write ^ Ai+ A2 with Ai = WI2WI1 and 
A2 = WliWl2- Apphcation of the rule in Eq. (EH) 
yields 

^1 = (\/Pl)2.2(V^)l,l' 

A2 = {Vp^),^,{Vp^)2,2 (40) 

and after using the rule in Eq. H34|l we obtain 
Ai = (VpT) 2,2(^/^)3,3 = v/a^, 

A2 = (vpr)^^^(ypr)i_i = Ai. (41) 

Thus, we obtain 

Hll = Ulul (Ai + Ta^) , (42) 

where we have used (— 1)™^"'^ det it| = t^i2f^2i- The re- 
maining V'f,l\ V^p], Vf^p^, Ppapi, T'pspa; ^p^pi are given 
by appropriate permutations of the A's. For to = 3 (full 
permutation) the only possible contributions are 

13(3) ^ 1 
P1P2P3 ' 

^^?)>3P2 = 3^/A^, (43) 



where we have used (—1)^ -'^detJ/" ~ +1 and, again 
from Eq. (|26|l . fpfp2P3 = 1 for A = 3. The expression 
for Vpip3p2 follows from Ai ~ A2 = A3 — VA1A2A3 and 
requires full rank to be nonvanishing. 

As an illustrative higher dimensional example we con- 
sider the case where A^ = 5, Z = 4 and a partial per- 
mutation specified by to = 2. The diagonal part can be 
calculated to 

^^tLp5P3 - t/|3^AiA3A4A5-KC/lUAiA2A4A5 

-KC/|5^AiA2A3A5 (44) 
and the permutation part is given by 

^^fUp3 = [(-l)'det4]V(fUp3. (45) 

where /p,p^p,p3 = Ai(Ai, . . . , A5) + A2(Ai, . . . , A5). For 

a calculation of fptlnp^pa we use the algorithm described 
above. From Eq. 1)30(1 we know that 

Ai = Wl2Wl,Wl2Wl,, 

A2 = Wl,Wl2Wl,Wl2. (46) 

Applying rule ((31|l to the indexes we obtain 

Al - (^)2.2(^)l.l(^)2,2(^)l,l' 

= (^)i.i(^)2.2(^)l,l(^)2,2- (47) 

After a transformation according to the rule H34|l we have 

Al = (^)2,2(^)L(^)m= VW^, 

= ^AiA4A2A5 (48) 
and consequently 

/pfp4P5P3 =Ai+A2^ v^A2AiA4 + ^AiA2A4A5. (49) 

Vpf'pipsps can now be written as 

^^'kp5P3 = (f^l24)' (v/A2AiA4+ ^AiA2A4A5) , 

(50) 

using det = — f/jj U21 ■ 

V. EXPERIMENTAL PROCEDURE 

When we consider possible experimental realizations 
of the off-diagonal mixed state phases we immediately 
encounter a problem: how do we experimentally imple- 
ment the Zth root of density operators? Fortunately, this 
may be resolved in the / = 2 case in the sense of purifi- 
cation, i.e., by adding an ancilla system in a certain way. 
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Here, we demonstrate this in general and propose a phys- 
ical scenario for the qubit case in terms of polarization- 
entangled two-photon interferometry. 

We first show how to realize the I — 1 and I — 2 phases 
via purification. For an N dimensional Hilbert space Ti., 
consider the nondegenerate density operator 



N 



Pi = ^>^k\lpk){lpk 



(51) 



fe=i 



A purification of this pi is any pure state ^'i obtained 
by adding an ancilla system a to the considered system 
s such that pi = TTa\'9i){'^i\. Thus, we may write 



l*i> 



N 

E 

k=l 



(52) 



where {jt/:'/!;)} is an orthonormal set of vectors in the an- 
cilla Hilbert space Ha- Consequently, any orthogonal 
density operator p„ 
tion of the form 



[Ug)" ^pi(C/J)" has a purifica- 



(53) 



for any unitarity Ua acting on Tia- In the following, we 
assume dimTia — N and put \ipk) — IV'fe)- 

Let Us (gi Ual^i) and Vs (g> Va\^i) be two Hilbert space 
representatives of a pair of purifications of UsPiUj and 
VsPiVg ■ The coincidence interference pattern obtained 
in superposition is determined by the interference profile 



Us(g>Ua\-^l) +Vs(g>Va\^l 



= 2 + 2Re[Tr(C/t^,®[/tv;|*i)(*i|)]. (54) 

By choosing = e'^{Ugy\ Vs = C/H (C/g)''"\ and 

Ua = Va — I, we obtain the / = 1 phase factors 7pJ] by 
variation of the U(l) phase x since 



^'^^[TT{{Uiy'u\Ugy^ < 



->^<i>[Tr(C/llp,J], 



where we have used that Tr^ [{Ug 



(55) 

->i>(*i|(t/t)^-^] 



(2) 

Pj^ . Similarly, the / — 2 phase factors ^p-^ p-^ are ob- 



tained by letting Us 

\ i2 - 1 



iUg) 



32-1 



Vs = U\^ (Ua 



Ua = (UgY' \ and Va = (C/ll)^(C/g)''' \ T being trans- 
pose with respect to the ancilla basis {iV'fc}}, since 

<^>[Tt{UIVs ®UlVa\^ l){^l\)] 

= e-'^<i>[Tr{{Uiy'u\\{Ugy' 

^{uiy'-\u^^f{Ugy^-'\^,){^,\)] 

= e-^'^$[Tr([/IIVp-[/IIVp-)], (56) 



where the last equality may be obtained by explicit use 
of l^-i) in Eq. (El with \Lpk) = IV'fe), V/c. 

Let us discuss a physical purification scenario for the 
I — 1 and / = 2 phases in the qubit case. Consider 
the two-photon Franson-type [2^ setup in Fig. |21 A 
source that in the horizontal-vertical {h — v) basis pro- 
duces polarization-entangled photon states of the form 



^{l + r)\h)®\h) 



+ \hi^-r)\v)®\v) 



(57) 



has been demonstrated in Ref. |2g. Considered as sub- 
systems both photons are in a mixed linear polarization 
state pi with polarization degree r. The desired super- 
position of Us (8) C/al^'i) and Vs (8) V^al^'i) is obtained by 
requiring sufficiently short coincidence window so that 
detection occurs only when the photons both either took 
the shorter path or the longer path |23|. A purification 
of the orthogonal density operator p2 = UgPiU^ may 
be achieved by flipping the polarizations of the photons, 
yielding 



\^2)=Ug®Ug\^,) = jhl+r)\v)®\v) 




+ \l^a-r)\h)(8>\h). (58) 



Us 



Vs 



FIG. 2: Franson setup for polarization-entangled photon 
pairs. In the longer arms, the system and ancilla photons 
are exposed to the polarization affecting unitarities Us and 
Ua, respectively, and similarly Vs and Va in the shorter arms. 

To demonstrate the I — 1 and I — 2 geometric phases 
in this scenario, it is sufficient to consider unitarities that 
rotate linear polarization states along great circles an an- 
gle (3 on the Poincare sphere, see Fig. This amounts 
to 



U{(3,e) = exp(-z^ 



cose{\h){v\ + \v){h\) 
z\h){v\+z\v){h\)]), (59) 



which fulfills the parallel transport conditions in Eqs. Q 
and Q with respect to the h — v basis. In practice, 
U (/3, 9) may be implemented by appropriate A— plates, 
the thickness and orientation of which correspond to the 
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parameters (3 and 0, respectively. For example, Ug = 
[/(tt, 7r/2) acting on the linear polarization states as a 
polarization flip and thus connects pi and p2 , is achieved 
by a A/2 plate with half axis making an angle 45° to the 
vertical (v) direction. Fm'thermore, 9 = and f3 = 7r/2, 
corresponding to a A/4 plate oriented along the vertical 
direction, takes h and v into the right (i?) and left (L) 
circular polarization states, respectively. 



h 




FIG. 3: Effect of the unitarity U (/3, 9) on the Poincare sphere. 
The horizontal polarization state h at the north pole is taken 
into a new polarization state at spherical polar angles {P,9). 

The phase factors 7pJ,\ n = 1,2, are obtained from 
the coincidence intensity by choosing Us = e'^^{Ug)^ , 
Vs = U{J3,e){Ugy^ , and Ua = Va = I. Explicit cal- 
culation for t/(/3,6i) in Eq. ^ yields Tt{piU {(3,9)) = 
Tv{p2U{(3,9)) = cos(/3/2), which entails that 7^^^ and 

7p2'' are real-valued and changes sign at /3 = (2j -f l)7r, 
j integer, corresponding to a sequence of phase jumps 
of TT. Furthermore, the choice Us = e'^^Ug, Vs = 

Ui/3,9), Ua = Ug, and K = U^i(3,9) = U{p,-9), 
(2) 

yields 7pip2 and we may compute the expected output 

as TT{^U{f3,9)^U{(3,9)) = x/T^cos^ (/3/2) - 

sin^ (/?/2), which is independent of 9 and can be posi- 

(2) 

five and negative for r 1 depending upon /?. 7pip2 
changes sign at/3 = 2Trj + 2 arctan \/l — r^. Note that 
< arctan v^l — < 7r/4, modulus tt, which assures 
that the I = 1 and I = 2 phases never become indetermi- 
nate for the same /3 value, and thus provide a complete 
experimental phase characterization of the qubit case in 
the sense of purification. 



VI. PROJECTION PHASE 

Since the definition of the off-diagonal geometric mixed 
state phase claims to be reducible to the pure state off- 
diagonal geometric phase, the question arises if there is 
a connection to the experimental verification of the lat- 
ter performed by Hasegawa et al. [TsL IT^ . This has 
to be answered in the negative, since in this experiment 
the evolution of the orthogonal state is implemented as a 
projection operator, which is by definition equivalent to 
a pure state. On the other hand, by taking the impurity 
of the input state into account, the shift in the interfer- 
ence pattern in the Hasegawa et al. setup, given by the 
additional phase factor 

7pP = $[Tr([/pC/P)], (60) 

could be used as a definition of the off-diagonal geomet- 
ric mixed state phase, if the unitarity U describing the 
evolution inside the interferometer is parallel transport- 
ing the eigenvectors of the nondegenerate p. Here, P 
should project onto the eigenstate that corresponds to 
the smallest eigenvalue of p. Parallel transport is for ex- 
ample fulfilled in the Hasegawa et al. experiment if the 
incident spinor is polarized in a plane perpendicular to 
the direction of the magnetic field. 7pp is gauge invariant 
under a U(l) transformation of each of the basis vectors 
and it reduces to the corresponding off-diagonal geomet- 
ric phase factor of Ref. ^3 in the pure state limit when 

u = u\\. 

In the two dimensional case relevant for the Hasegawa 
et al. experiment with input p = Ai|-!/'i)(V'i| + '^2|'02)(V'2|j 
Al > A2, we can write Eq. as 

TrlUpUP] = Ai(-l + 772) + A277'e"2'". (61) 

Here, U G SU(2) with the diagonal matrix elements 
Uii = U22 = '76*" is not necessarily fulfilling the par- 
allel transport condition with respect to |^/'2)}. In 
the pure state limit Ai = 1, A2 = the off-diagonal 
phase is always tt since Tr[UpUP]\^^i = — 1 + ry^ is real 
and negative, irrespective of whether U parallel trans- 
ports \4'i)t\^'2) or not. For a mixed input state p the 
A2-term does not vanish and we obtain additional ge- 
ometric and/or dynamical phase contributions. These 
can be considered to originate in the subjacent geometry 
only if J7 is a parallel transporting unitarity, but not for 
arbitrary U. 

To show the consistency with the experiment per- 
formed by Hasegawa et al. we calculate the phase 
(f>pP = arg Tr [C/p[/P]. In the left panel of Fig. ^we show 
(j)pP for a mixed input state with Ai = 0.87, A2 = 0.13, in 
accordance with the experimental degree of polarization 
in Refs. 15, 16], and the spin polarization angle 9 = tt/6 
relative to the magnetic field in the upper arm of the 
interferometer (see Fig. 2 of Ref. |0|). The calculated 
curve matches with the experimental and theoretical re- 
sults presented in Fig. 5(d) of Ref. Note that in this 
case U is not parallel transporting the incident spinor. 
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Another interesting fact is that due to the impurity of 
the input state we expect phase jumps for 9 — -k /2 for 
(5 = 2 arccos and 5 — 2t: — 2 arccos see right 
panel of Fig. ^ where 5 is the precession angle of the 
incident spinor about the direction of the magnetic field. 
Here, we have a parallel transporting U = C/", thus these 
jumps have their origin in the subjacent geometry of state 
space. 
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FIG. 4; The projection off-diagonal mixed state geometric 
phase arg Tr[C/pC/P] in radians modulus 27r for the Hasegawa 
setup with Ai = 0.87, A2 = 0.13 and Q = 7r/6 or = 7r/2, re- 
spectively . 5 is the precession angle in radians of the incident 
spinor about the direction of the magnetic field. 

The projection off-diagonal geometric phase factor 7pp 
is invariant under phase transformations of the eigenvec- 
tors of /o, it reduces to the corresponding I ~ 2 phase 
factor of Ref. [l3 in the pure state limit, and it has the 
advantage that it can be observed in single particle exper- 
iments. The drawbacks are that it is less symmetric than 
that in Eq. one cannot easily state a generalization 

like in Eq. (|14|l . and vice versa it cannot be regarded an 
off-diagonal generalization of the mixed state geometric 
phase in Ref. P| as p and P are not unitarily connected. 
These features suggest that the mixed state geometric 
phase factor in Eq. (|16|l is to prefer over 7pp. 

VII. CONCLUSIONS 

Recent investigations in geometric phases in quantum 
systems have led to cases where the standard definitions 



breaks down. On one hand, such situations emerge for or- 
thogonal initial and final pure states connected unitarily, 
on the other, unitary evolution of a system in a mixed 
state may lead to nodal points in parameter space. In 
search for a complementary geometric quantity defined in 
such cases, the off- diag onal mixed state geometric phase 
has been proposed (iJl by the present authors as a gen- 
eralization of the off-diagonal geometric phases for pure 
state put forward in Ref. ^3 ■ 

Starting with a preliminary discussion about orthog- 
onality of mixed states and quantum parallel transport 
we have provided a general treatment of the off-diagonal 
mixed state geometric phase comprising unitarities that 
can be decomposed into a diagonal part leaving the initial 
basis states unchanged and a permutation part reorder- 
ing the initial states. An algorithm has been presented 
to calculate the appropriate phase factors efficiently for 
any dimension and for an arbitrary number of orthogo- 
nal density operators. Furthermore, we have discussed 
the projection off-diagonal geometric phase appearing in 
the neutron experiment by Hasegawa tt al. ,15, 16] as an 
alternative definition of off-diagonal mixed state phase. 

In the qubit case the off-diagonal mixed state phase 
can be fully qualified both from the theoretical and from 
the experimental point of view. But it has to be men- 
tioned that the measurement seems to require control and 
measurement of one or more ancilla systems although the 
off-diagonal mixed state phases are properties of the sys- 
tem alone, since the constituting set of density operators 
pertains solely to the system. Explicitly, a Franson in- 
terferometer setup for the qubit case has been presented 
illustrating the nontrivial sign change property of the off- 
diagonal phase connected to the mixed state case. The 
apparent need for control over an ancilla system seems to 
suggest that the proposed concept of off-diagonal mixed 
state geometric phase is a nonlocal and/or contextual 
property of the unitary evolution of a quantum system. 
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